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Bottomonium spectroscopy with mixing of rjb states and a light CP-odd Higgs 
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The mass of the r]b(lS), measured recently by BABAR, is significantly lower than expected from 
QCD predictions for the T(IS') - rjt(lS) hyperfine splitting. We suggest that the observed r)i,(lS) 
mass is shifted downwards due to a mixing with a CP-odd Higgs scalar A with a mass itia in the 
range 9.4 - 10.5 GeV compatible with LEP, CLEO and BABAR constraints. We determine the 
resulting predictions for the spectrum of the rjt(nS) — A system and the branching ratios into r + r~ 
as functions of ntA- 
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The BABAR collaboration has recently determined the 
rjb(lS) mass m^ns) with an error of only a few MeV in 
radiative decays T — > jrfo of excited T states and the 
observation of peaks in the photon energy spectrum. The 
result from T(35) decays is m^is) — 9388. (stat)± 
2.7 (syst) MeV 0, and the result from Y (25) decays is 
m Vb(1S ) = 9392.9±t;| (stat) ± 1.9 (syst) MeV 0. The 
average gives 0] 



"^6 (is) 



= 9390.9 ±3.1 MeV 



(1) 



implying a hyperfine splitting Ehf s (lS) = rnyns) 
"V(is) of 



Elf s (lS) = 69.9 ±3.1 MeV 



(2) 



This result can be compared to predictions from QCD. 
Recent results based on perturbative QCD are in good 
agreement with each other and give Ehf s (lS) = 44 ± 
11 MeV and E hfs (lS) = 39 ± 14 MeV Q (whereas 
Ehfs (I'SjL, varies over a wider range in phenomenological 
models Q); in the following we will use tentatively an 
average value 



E hfs D ( lS ) = 42 ± 13 MeV 



(3) 



which is about two standard deviations away from the 
experimental result ([2]). The most recent result from (un- 
quenched) lattice QCD is [f| 



E 



latQCD 
hfs 



(15) = 61 ± 14 MeV 



(4) 



which is within 1 a of ©. However, the hyperfine split- 
ting is quite sensitive to short distances or hard quark 
momenta. It has been argued in Q that the perturbative 
results in Q can be used for short distance corrections 
of the Wilson coefficient of the corresponding spin-flip 
operator measured on the lattice. The additional con- 
tribution 5 hard E hfs (lS) to E l ^ff CD (lS) has been esti- 
mated as S hard E hfs (lS) 20 MeV 7], which brings 

the lattice result dl} in good agreement with the pertur- 
bative result ([3]). Although this conclusion needs to be 
checked within perturbation theory with lattice regular- 
ization, we consider it as a support for the perturbative 
QCD result ©. 

Whereas an explanation of the discrepancy between 
@ and © within QCD is not excluded at present, we 



will elaborate below the consequences of an explanation 
of this discrepancy due to new physics in the form of a 
mixing of the % states with a CP-odd Higgs scalar A with 
a mass itia in the range 9.4— 10.5 GeV: as a result of such 
a mixing, the masses of the rjb - like eigenstates of the full 
mass matrix can differ considerably from their values in 
pure QCD without the presence of the CP-odd Higgs A 
H B El , and the mass of the state interpreted as the 
r)b (IS) can be smaller than expected if wia is somewhat 
above 9.4 GeV. 

In such a scenario, the masses of the states interpreted 
as r]b(2S) and r]b(3S) can also be affected, and all states 
can have non-negligible branching ratios into t + t~ due 
to their mixing with A. According to recent results of 
BABAR [11(, the corresponding branching ratio of the 
observed state is below 8% at 90% confidence level. The 
branching ratio into fi + fjT would be smaller by a factor 
m^/m 2 , and well below the present upper limit [l2|. The 
investigation of these phenomena is the purpose of the 
present paper. 

A relatively light CP-odd Higgs scalar can appear, e. g., 
in non-minimal supersymmetric extensions of the Stan- 
dard Model (SM) as the NMSSM(Next-to-Minimal Su- 
persymmetric Standard Model) Cfl, El El ES EE IS • 
Its mass has to satisfy constraints from LEP, where it 
could have been produced in e + e~ — * Z* — * Z H and 
H — > A A (where H is a CP-even Higgs scalar). For 
itia > 10.5 GeV - where A would decay dominantly into 
bb - and uih < 110 GeV, corresponding LEP constraints 
are quite strong [l8j . For 2m T < rriA < 10.5 GeV, A 
would decay dominantly into t + t~ and values for mjj 
down to ~ 86 GeV are allowed [l8| even if H couples to 
the Z boson with the strength of a SM Higgs boson. 

In fact, searches for e + e~ — > Z* — > Z H with H — > bb 
indicate a light excess of events (of ~ 2.3 a significance) 
for ma ~ 95 — 100 GeV [l8j], which could be explained 
by a reduced branching ratio BR(H — > bb) ~ 0.1 
and a dominant branching ratio BR(H -t 44) ~ 0.9 
EE E3 if A decays dominantly into r + r~ . The pos- 
sible explanation of this excess of events at LEP is an 
additional motivation for a CP-odd Higgs scalar with a 
mass below 10.5 GeV. Allowing for m# somewhat below 
100 GeV, such a scenario would also alleviate the "little 
fine tuni ng pro blem" of supersymmetric extensions of the 

SM [HEllil- 



2 



Finally CLEO and BABAR [n| have searched for 
a light CP-odd scalar ^4 with A — > r + r~ in radiative 
T(15) and T(35) decays respectively, which mainly con- 
strains the range ttia < 9.4 GeV. In this work we focus 
on the range 9.4 GeV < rriA < 10.5 GeV, which would 
be the most relevant for strong ly, — A mixing effects as 
advocated in (Toj . 

Apart from ttia, such mixing effects depend in a 
calculable way on the model-dependent coupling of A 
to b quarks. The corresponding coupling, normalized 
with respect to the coupling of the SM Higgs scalar 
to b quarks, will be denoted by X^. In models with 
two Higgs doublets H u and Hd, where H u couples to 
up quarks and Hd to down quarks and leptons, one has 



Xa — cos ( 



tan j3 



(5) 



where cos 9a denotes the SU(2) doublet component of the 
CP-odd scalar A, and tan/3 = (H u ) / (H d ). For tan/3 > 
1, Xd can equally satisfy Xd 1 |10| . 

Below we proceed as follows: i) First we assume that, 
in the absence of a CP-odd Higgs scalar, m^o^gj would 
have a value compatible with ([3]), 



i. e. 



9418 ± 13 MeV 



(6) 



ii) We diagonalise the rjb(nS) — A mass matrix (n — 
1, 2, 3) and require that one eigenvalue coincides with the 
mass measured by BABAR within errors ([T]); this con- 
dition gives us an allowed strip in the Xd — vtla plane, 
which depends only weakly on the assumed masses of 
7^(25) and 7,° (35). 

The resulting values of Xd as function of ttia allow to 
determine the remaining eigenvalues of the mass matrix, 
and the decompositions of the eigenvectors in terms of A 
and r]®(nS). Finally the A components of the eigenstates 
allow us to determine their partial widths and to estimate 
the branching ratios into t + t~ . 



For the 7^(15) 
make the ansatz 



r)®(2S) — T]®(3S) — A mass matrix we 



M 2 = 



( 



\ 






5m\ 





Sm?, 






%°(3S) 

Sm 2 



8m\ \ 



5m 2 
Sm 2 



(7) 



/ 



The diagonal elements m 2 a , „, are assumed to be known 
(within errors) from QCD with ^ for m^ons)- The 
results depend only weakly on the ^(25) and rj® (35) 
masses, for which we take [|[ m^o^s) = 10002 MeV, 
77^^0(35) = 10343 MeV. (We neglect finite width effects in 
A4 2 , and take real matrix elements.) 

The off-diagonal elements bm 2 n can be computed in 
the framework of a non-relativistic quark potential model 
in terms of the radial wave functions at the origin 0, 
9]. These can be considered as identical for vector and 
pscudocalar states, and be determined from the measured 
T — ► e + e~ decay widths. Substituting recent values for 



these widths (see [l(| for details) we find 
5m\ 



8m\ 



(0.14 ± 10%) GeV 2 x X d , 
(0.11 ± 10%) GeV 2 x X d , 
(0.10 ± 10%) GeV 2 x X d . 



(8) 



We estimated the errors from higher order QCD correc- 
tions to the relation between the radial wave functions 
at the origin and T(T(nS) -> e+e") to be - 10%. 

(These errors play only a minor role for our results.) 

In order that one eigenvalue of M 2 coincides with the 
BABAR result ([1]) subsequently denoted as m b s , niA in 
([7]) has to satisfy 



8m\ 



obs 



in 



r,°(XS) m obs 



8m\ 



<5to| 



X(2S) 



S«(3S) 



(9) 



obs 



Once m,A is expressed in terms of X d , X d remains the 
only unknown parameter in M 2 . Varying m f, s within 
the errors in ([T]), m^o^g^ within the errors in ^ and Sm 2 
within the errors in ([5]), we obtain for Xd as a function of 
iriA the result shown in Fig.[T] (Recall that wla - or, the 
heaviest eigenstate of M 2 with a large A component - has 
to be below 10.5 GeV in order to satisfy LEP constraints 
in the presence of a Higgs scalar with a mass around 
100 GeV.) 




m. (GeV) 1 



FIG. 1: Xd as a function of ttia (in GeV) such that one 
eigenvalue of .M 2 coincides with the BABAR result Q. 

Now the masses of all 4 eigenstates of A4 2 can be com- 
puted, which are shown together with the error bands 
from m bs, TO >jO(i5) an< i (i n orange/grey) in Fig. [5] 
as functions of ttia- Henceforth we denote the 4 eigen- 



states of A4 by r]i 



1 ... 4 where, by construction, 



m b s ■ For clarity we have indicated in Fig. [2] our 
assumed values for m^o^m as horizontal dashed lines. 
For m,A not far above 9.4 GeV (where Xd is relatively 
small) the effects of the mixing on the states 77° (25) and 
77°(35) are negligible, but for larger itia the spectrum 
can differ considerably from the one expected without 
the presence of A. 

Now we consider the branching ratios of the eigenstates 
into t + t~ , which are induced by their A-components. 
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FIG. 2: The masses of all eigenstates as function of tua- 



The decomposition of the eigenstates into the states be- 
fore mixing can be written as 

m = P M rfi(lS) + P it 2 %°(25) + P i>3 r? 6 °(35) + P M A . 

(10) 

It turns out that the coefficients Py can be expressed 
analytically in terms of the eigenvalues m 2 . of M 2 : 



+ 



5mk 



8m\ 



Pi 



—5m 



1,3 



-P, 



i,4 



for j = 1,2, 3 



-1/2 



(11) 



In Fig. [3] we show our results for the A-components 
Pi. 4 for all 4 eigenstates together with the error bands 
from m i, s , m J) °(is) and 5m 2 . 

In the case of rji = rj obs , only the coefficients Pi.i and 
Pi ,4 differ significantly from 0. This allows to express the 
branching ratio BR(r/i 



t + t ) as 



BR( m 



-) 



p2 rrT 
^M 1 A 



where T A T is the partial width for A 



the width of the state r/ b (lS) without mixing with A, and 
T^ 4 the total width of A (without mixing). For T A T we 
have 



-pTT 

1 A 



m 



4 JP (13 > 



where we would have to identify M with if A decays 
on its mass shell. Since the dependence on M originates 
from phase space integrals, M has to be identified with 
the mass of the decaying state (which is actually always 
close to toa) in our case. Since the BR(A — > r + r~) is 
typically - 90% (for large tan/3), we take r*|* ~ 1.1 T T A T 
in Cr 




m 4 (GeV) 



FIG. 3: The A-components |Pi,4| for all 4 eigenstates as func- 
tions of 171 a • 



It is remarkable that, once we insert eqs. (fTTj) and (fl3|) 
into the expression (fT2")l for BR(rji — > t + t~), all depen- 
dence on cancels: even for — » 9.4 GeV where 
Pf 4 — » 1, one has r^* — > from Xd — » 0, and hence the 
first term in the denominator in (|12p always dominates. 
The numerical result is 



BR{ m 



') = {2A+_ 2 2 D x 10- 2 x 



lOMeV 



.(14) 

Since the analytic expression for this BR is proportional 



to (m 2 ob 



m 



v ms)> 



in our approach, its lowest possible 



value is quite sensitive to the smallest allowed value for 
this difference. 

Using T v o(nS)/T Va (nS) ~ (m b /m c )[a s (m b )/a s (m c )} 5 
~ 0.25 - 0.75 (2i and T rh (lS) = 26.7 ± 3 MeV [p 



we estimate T 



»?? (IS) 



5 — 20 MeV. Hence the predicted 



branching ratio is compatible with BABAR upper limit 

of8%[nf. 

Turning to the remaining heavier eigenstates, expres- 
sions similar to (| 12[) are always very good approximations 
for the branching ratios into r + t~ , since the eigenstates 
consist essentially of just one rj b (nS) state and the CP- 
odd Higgs A. However, the branching ratios vary now 
with m rji and hence with mA', the results are shown 
(12) graphically in Fig. 2] assuming r^o^g) ~ 10 MeV and 
r»)0(2S) ~ r»;{J(3S) ~ 5 MeV. For larger (smaller) total 
widths these branching ratios would be somewhat smaller 
(larger) . 

An important issue is the production rate of the eigen- 
states rji in radiative decays of excited T states, notably 
in T(3S) — > 777^. For a pure CP-odd scalar A, the 
BR(T(3S) -> 7 A) is given by the Wilczek formula [H 



BR{T(nS) ^-fA) 
BR (T(nS) 



ji+ji ) 



\/2vra 



l T(3S) 



x F 



(15) 



F is an dependent correction factor, which includes 
three kinds of corrections (the relevant formulas are sum- 
marized in [23j): bound state, QCD, and relativistic cor- 
rections. Unfortunately these corrections become unreli- 
able for mA > 8 GeV; a naive extrapolation of the known 
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FIG. 4: The branching ratios into r r for the eigenstates 
r/2, r/3 and r/4 as functions of uia. 

corrections leads to a vanishing correction factor F for 
rriA > 8.8 GeV as relevant here. 

Thus it is difficult to predict the branching ratios 
BR(T(3S) -► 7^): if the BR (T(3S) -► 7 A) is assumed 
to vanish, the production of the states has to rely on 
their r)®(nS) components. Otherwise, negative interfer- 
ence effects could appear leading to suppressed branching 
ratios. Hence it cannot be guaranteed that the (kinemat- 
ically accessible) part of the spectrum shown in Fig. [2] is 
actually visible in radiative T(35) decays in the form of 
a peak in the photon energy spectrum . 

In view of the possibly quite low photon energies 
and/or large backgrounds, the photons can well escape 
undetected even if the process T(35) — ► 7^ occurs with 



a non-negligible rate. In this advocated in 24], 

the A components of rji can still manifest themselves 
through a breakdown of lepton universality in the form of 
an excess of t + t~ final states in T(35) — * l + l~ fTol] . 
However, in the case of a vanishing correction factor F 
in (|15[) . this phenomenon would disappear as well. 

On the other hand, if the existence of a CP-odd scalar 
A and a CP-even scalar H with mu ~ 95 — 100 GeV 
(and a dominant H^AA branching ratio of ~ 90%) 
is responsible for the excess of events at LEP as noted 
above, it becomes important to test this scenario at the 
LHC: The standard search channels for a SM-like CP- 
even scalar H would fail, and the final states from H — > 
A A with mA below the bb threshold would be difficult 
to detect. Proposals for a verification of this scenario at 
the LHC have been made recently in [2^, [2(| [53] • 

To conclude, if the mixing with a CP-odd Higgs scalar 
is responsible for the discrepancy between the BABAR 
measurement of m Vb (is) and the expectations from QCD, 
it can manifest itself in the form of a completely distorted 
spectrum of states as shown in Fig. [51 The branching 
ratios into t + t~ would be non- vanishing, albeit below 
the present experimental upper limit (for the lowest lying 
state). These manifestations of a light CP-odd Higgs 
scalar in T physics at (Super) B factories [28[ would be 
complementary to its possible discovery at the LHC. 
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